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ABSTRACT 



One prediction of particle acceleration in the supernova remnants in the magnetic 
wind of exploding Wolf Rayet and Red Super Giant stars is that the final spectrum is a 
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composition of a spectrum E~ 7 I 3 and a polar cap component of E~ 2 at the source. This 
polar cap component contributes to the total energy content with only a few percent, but 
dominates the spectrum at higher energy. The sum of both components gives spectra 
which curve upwards. The upturn was predicted to occur always at the same rigidity. 
An additional component of cosmic rays from acceleration by supernovae exploding into 
the Inter-Stellar Medium (ISM) adds another component for Hydrogen and for Helium. 
After transport the predicted spectra J(E) for the wind-SN cosmic rays are E~ 8 ^ 3 and 
E~ 7 / 3 ; the sum leads to an upturn from the steeper spectrum. An upturn has now been 
seen. Here, we test the observations against the predictions, and show that the observed 
properties are consistent with the predictions. Hydrogen can be shown to also have a 
noticeable wind-SN-component. The observation of the upturn in the heavy element 
spectra being compatible with the same rigidity for all heavy elements supports the 
magneto-rotational mechanism for these supernovae. This interpretation predicts the 
observed upturn to continue to curve upwards and approach the E~ 7 / 3 spectrum. If 
confirmed, this would strengthen the case that supernovae of very massive stars with 
magnetic winds are important sources of Galactic cosmic rays. 

1. Introduction 

Since 2008 there has been increasing evidence for an extra component of cosmic ray electrons 
and positrons, with a plateau of cosmic ray electrons showing an E~ 3 spectrum, and the CR- 
positron/electron ratio approaching E +1 ^ 3 (ATIC: Chang et al. 2008; Pamela: Adriani et al. 
2009). These results confirm a quantitative model originally proposed in 1993 (Stanev et al. 1993), 
in which the supernova shock racing through a magnetic wind (Biermann 1951, Parker 1958) gives 
rise to a source spectrum of E~ 2 from the polar cap of the star, where the magnetic field is radial, 
and E~ 7 1' 3 from most of 4ir, where the magnetic field is nearly circular (Biermann 1993). Allowing 
for losses (Kardashev 1962) the polar cap component spectrum becomes E~ 3 for observers of cosmic 
ray electrons. Since acceleration is slower for a radial magnetic field (Jokipii 1987; Meli & Quenby 
2003a, b; Ellison & Double 2004; Meli & Biermann 2006) than for a tangential magnetic field, more 
secondaries are produced in the polar cap region, any resulting secondary spectrum is shifted down 
in energy and up in flux thus explaining the observations (Biermann et al. 2009b). 
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The prediction of the E~ 3 cosmic ray electron component was confirmed (HESS: Aharonian et 
al. 2009) with a spectrum of E~3.i±o.i(stat.)±OA(syst.) _ rjij^ g a pp roa ch has also led to an understand- 
ing of the WMAP haze (Finkbeiner 2004; Hooper et al. 2007) with a predicted spectral behavior 
and v 3 / 2 (Biermann et al. 2010). This prediction is consistent with the detection 
of the Inverse Compton component of the WMAP haze, the Fermi haze and its spectrum (Dobler 
et al. 2010; Su et al., 2010). These spectral index predictions are derived from the assumption, 
that a diffusive regime in a disk connects to a convective regime in a magnetic wind for cosmic 
ray transport. Thus, several independent observations support the quantitative model proposed in 
1993 (Stanev et al. 1993) of massive star explosions. The original cosmic ray spectra are consistent 
will all these new observations. 

There are a number of ideas that have suggested an upturn in the cosmic ray spectra, such as 
a nearby source or substructure in the accelerating shocks. Here, we focus on the upturn of the 
spectra predicted in 1993 (Biermann 1993, Stanev et al. 1993, parallel papers are summarized in 
Biermann 1994). 

It had been proposed (see also Prantzos 1984, 1991, 2010; Vdlk & Biermann 1988) that the 
origin of cosmic rays can be traced to three source components: a) The cosmic rays originating 
from supernovae exploding into the ISM; b) The cosmic rays from supernovae exploding into red 
supergiant (RSG) winds; and c) The cosmic rays from supernovae exploding into blue supergiant 
or Wolf Rayet (WR) star winds. Disregarding the possible binary character of the stars, this is 
essentially a mass sequence, and a key property of these stellar winds is that they are magnetic 
(Abbott et al. 1984; Barvainis et al. 1987; Drake et al. 1987; Churchwell et al. 1992). As Parker 
(1958) showed, the asymptotic magnetic field topology is radial in a polar cap, and tangential over 
most of 4tt. Furthermore, observations of these stars and their interpretation (e.g. Woosley et al. 
2002, Heger et al. 2003) show that the wind eats back down into the star, and so the abundances in 
the RSG star winds are enriched in Helium, and those of WR stars are concentrated in Carbon and 
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Oxygen. Therefore it can be expected, that in cosmic rays at higher energy much of the Helium 
can be traced to these winds of RSG stars, and possibly all of the elements heavier than Carbon 
to the WR stars (Stanev et al. 1993, Biermann 1994). 

We use here the concept that interstellar transport (e.g., Rickett 1977) is governed by a Kol- 
mogorov spectrum (Kolmogorov 1941a, b, c), as discussed in Biermann (1993), and Biermann et 
al. (2001). We test this prediction with new data (CREAM: Seo et al. 2008; Ahn et al. 2009, 
2010). The data are shown to be consistent with the prediction of 1993, and so we propose further 
tests. 

2. The upturn observed in cosmic ray spectra 

The general concept developed in 1993 is that, after transport is allowed for, the spectra are 
composed of a sum of a E~ s ^ 3 component, and a E~ 7 / 3 polar cap component, which is at a level of 
a few percent in integrated power. These components run up to the knee energy, with about Z 10 15 
eV, where the polar cap component cuts off completely, and the main component turns steeper, at 
about E~ 3,15 (Biermann 1994, Biermann et al. 2005). The spectrum finally cuts off completely at 
about Z 10 17 eV. Here Z is the charge of the cosmic ray nucleus. 

It is not obvious a priori, whether RSG star winds should have the same properties in terms 
of their magnetic field as WR star winds: It is not clear, as to whether the energy fraction of the 
polar cap component is the same for both, and also whether the characteristic knee rigidity ratio is 
the same. However, for simplicity we will assume that the population energy fraction and the knee 
particle energy are the same for both kinds of stars until the data force us to accept otherwise. 

We first deal with all elements from Carbon (C) through Iron (Fe) and with the combined data 
as a function of energy/nucleon. Next we deal with Helium, and test whether the same numbers 
can describe its spectrum. Finally we deal with Hydrogen, which does have a major ISM-SN 
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contribution. In our 1993 prediction Hydrogen had a cross-over with Helium, when the spectra are 
plotted as a function of energy per particle, which is confirmed now. We also test whether Helium 
or any other elements also require an ISM component to understand its spectrum. 

2.1. Carbon through Iron 

5 /2 

The Ahn et al. (2010) data are plotted as E n J(E n ), where E n is the energy per nucleon. 
Since the prediction is that J(E n ) is the sum of one component with E n S ^ 3 and the polar cap 
component with E n 7 ^ 3 , we note that the product E^J 2 (E n 8 ^ 3 + E n 7 ^ 3 ) is a symmetric shallow 
near-parabola, with the break at the energy E n ^ m , where the two contributions have equal flux. 
Fitting the shape of this function involves one parameter, that energy per nucleon at the break 
energy E n>m . 

First of all, the spectra before the break are determined by CREAM (Ahn et al. 2009) to 
be £;-2.66±o.04^ to be compared with the 1993 prediction of ^;-8/3-o.02±o.02 , with an asymmetric 
predicted error distribution. The 2009 determination is reasonably close to the prediction, as the 

—8/3 

difference is only 0.007 plus the corresponding errors. So we stick with E n . The cosmic ray 
electron data suggest in turn, that at source E~ 2 is a good fit for the polar cap component, and so 

— 7/3 

we also stick with E n after transport. 

Fitting this to the curves of the elements reveals that the break energy E n ^ m is consistent with 
being the same for all these elements. 

The fit curve En" 2 (E n 8 ^ 3 + E n 7 ^ 3 ) to the plot demonstrates, that a break energy of 10 TeV/N 
is consistent with the data. The present data are in fact compatible also with 3 or 30 TeV/N to 
within the errors, but 10 TeV/N is the best common fit. The 1993 arguments gave 12 TeV/N, 
derived from a fit to overall air shower data, at vertical and slanted zenith angles. 

We show in Fig. Q] the fits of all elements heavier than Hydrogen compared to experimental 
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data. 

2.2. Helium and Hydrogen 

Measured in flux per energy/particle Helium crosses over Hydrogen near about 10 TeV. The 
specific spectral index of the ISM-contribution has to be determined best by the Hydrogen com- 
ponent, and that approaches ^;- 2 - 78 ± - 009 for protons at low energy (AMS: Alcaraz et al. 2000). 
Since there is a slight flattening due to the wind-component, we need to adopt a spectrum for the 
ISM-component, originally predicted to be E~ 2 ' 75:iz0 ' 04: , for which we write more generally E~ PlSM . 

Fitting then En" 2 (E n 8 ^ 3 + E n 7 ^ 3 + bE n PlSM ) gives us an estimate of the ISM contribution. 
We show the fit for H in Fig 12.21 and for comparison repeat the fit for He. 

Fitting the data (from AMS: Alcaraz et al. 2000) gives at higher energy a spectral index of 
2.66 ± 0.02, This suggests that the ISM contribution has pisM — 2.80 and then corresponds to 
bs — 1 for Hydrogen at the steeper spectral index 2.80, and about bjj — 4 for the slightly flatter 
index of 2.75. This corresponds to a wind-contribution of the energetics for CR-Hydrogen of 25 
% for the steeper index, and about 12 % for the slightly flatter index. Using pism = 2.70 does 
not give a satisfactory fit to the Hydrogen data. The differences illustrate the uncertainties. For 
Helium the ISM-CR component is very much smaller, and hard to define within the errors, using 
data from independent instruments. 

The largest modification with respect to 1993 is the slightly larger contribution of the wind 
component to Hydrogen. 

2.3. Differential spectra 

Differential spectra are shown in Ahn et al. (2009). All of those show near constant ratios, 
with the exception of Nitrogen: The Nitrogen/Oxygen ratio shows a weak dependence, consistent 
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Fig. 1. — Fit for elements heavier than H with the spectral shape from Stanev et al. (1993). We use 
for the break energy 10 TeV/N and only fit the fluxes above 50 GeV/N because of solar modulation 
and low energy spallation effects. No ISM-SN CR-component is necessary to fit the spectra. Lower 
energy data implies a small ISM-SN CR component for He. 
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Fig. 2. — Fit to proton and Helium data with spectral shape given in Stanev et al. (1993). For He 
we use for the break energy 10 TeV/N, and for H 20 TeV, since the rigidity E/Z is given. For H the 
strongest component is from ISM-SN, a wind-component is required to fit the data, but is weaker. 
Replotting these spectra in energy per particle would shift He up by a factor of 8 and to the right 
by a factor of 4, at 10 TeV/N, and so give the cross-over near that energy with Hydrogen. 



- 9 - 

with i?~ 5 / 9 suggesting interaction in an environment of a massive wind-shell, where the spectrum 
of magnetic irregularities excited by cosmic ray particles themselves dominates over other spectra 
(Biermann 1998; Biermann et al. 2001, 2009b). 

The Iron/Oxygen ratio shows a weak increase consistent with weak differential spallation: Low 
energy Fe has stronger spallation, so less iron is left. If this is the true interpretation, then at higher 
energy the ratio should approach a constant. The numbers suggest that up to about half of the 
low energy iron is lost through spallation. 

2.4. Supernova mechanism 

The rather remarkable observation, that this break energy E nm , where the E~ 8 / 3 and the 
polar cap component E~ 7 / 3 have equal fluxes, is well defined, even if not precisely known, imme- 
diately shows that all participating supernovae appear to approach asymptotically a commonality 
of correlated magnetic field strength and explosion energy (Biermann 1993). Stellar evolution of 
massive stars (see, e.g., Woosley et al. 2002) suggests, that indeed a commonality of properties 
in their asymptotic state might well be quite plausible. Using the Parker limit (Parker 1958) for 
a tangential azimuthal magnetic field this requires B^r) rUg N to be the same for all such su- 
pernovae, and since in this limit B<j,{r)r is a constant, using the base of that regime r s requires 
-S(/>,s fs U$ N to be the same; here Usn is the shock speed of the supernova racing through the wind, 
and ~ 1/r is the azimuthal magnetic field. Assuming that the mass ejected is similar this con- 
nects the explosion energy directly to the magnetic field strength. The only mechanism that does 
that is the magneto-rotational mechanism (Kardashev 1964, Bisnovatyi-Kogan 1970, Bisnovatyi- 
Kogan & Moiseenko 2008, Moiseenko et al. 2010). One could speculate that the neutrinos couple 
to the exploded magnetized outflow through instabilities and irregularities, and so enhance the 
explosion just as photons can couple better to a wind in early type stars since the wave speed in a 
magnetic wind is higher (Seemann & Biermann 1997); thus it is conceivable that neutrinos play a 
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concomitant role in the explosion mechanism (see, e.g., Woosley et al. 2002). 

Since this mechanism employs rotation and magnetic field to their limits, this mechanism also 
gives rise to the axial symmetry seen in some supernova explosions and required by Gamma Ray 
Bursts (GRBs). So it is tempting to speculate that this is in fact the same mechanism that leads 
to supernovae and GRBs. 

2.5. Higher cosmic ray energies 

The sum of all these spectra should obey the observed total spectra of particles determined 
in other ways. This should be consistent with the overall spectra all through EeV energies. This 
test was used in 1993, and now, with increased H and He wind-contributions, can be done again. 
Using the new data increased the Hydrogen and Helium wind contribution, and so we may have to 
decrease the break energy by a small amount to compensate in the total from the knee and beyond. 

We note that the particles near the knee, here with a notable contribution also from Hydrogen 
and Helium, may constitute the seed particles for the next step up in energy, when a relativistic 
shock (Gallant & Achterberg 1999) may push the entire spectrum up by a considerable factor. 
Ultra high energy cosmic rays could be drawn from a strongly enriched composition (Biermann 
et al. 2009a, Gopal-Krishna et al. 2010). With such a concept the results here suggest that in 
different parts of the sky, ultra high energy cosmic rays might be dominated by low mass and very 
high mass elements (Abbasi et al. 2009; Abraham et al. 2010). 

3. Conclusions 

Here we demonstrate that the upturn recently observed in cosmic ray spectra matches the 
quantitative predictions from 1993. As now a number of independent observations (CR-positrons, 
CR-electrons, WMAP haze, Fermi haze, the 511 keV line emission from the Galactic Center region, 
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and the CR-upturn observed by CREAM) have been shown to be consistent with the original 
quantitative proposal (Stanev et al. 1993), the implications originally suggested find support as 
well. 

Since the upturn is consistent with being at the same energy/charge ratio for all heavy elements, 
and clearly very nearly the same for the very large number of participating supernovae, this finding 
supports the magneto-rotational mechanism for massive star explosions (originally proposed by 
Bisnovatyi-Kogan 1970; see also Biermann 1993). Only if energy and magnetic field of the exploding 
star are strongly correlated, can such a characteristic feature in the observed cosmic ray spectrum 
be the same for all such stars, whose explosions contribute to the cosmic rays at high energy. 

A clear prediction is that the upturn should continue to approach asymptotically a spectrum 
of E~ 7 I\ 

As so often, much better data, extending to yet higher energy, will allow to refine and test the 
concept at yet higher accuracy. 

4. Acknowledgements 

PLB would like to thank G. Bisnovatyi-Kogan, L. Clavelli, T.K. Gaisser, B. Harms, A. Heger, 
P. Joshi, K.H. Kampert, Gopal-Krishna, N. Langer, S. Moiseenko, B. Nath, G. Pavala§, B. Sadoulet, 
E. Salpeter, R. Sina, V. de Souza, P. Wiita, and many others for discussion of these topics, and O. 
Ofoha for help of the analysis. Support for work with PLB has come from the AUGER membership 
and theory grant 05 CU 5PD 1/2 via DESY/BMBF and VIHKOS. JKB and JD acknowledge the 
support from the Research Department of Plasmas with Complex Interactions (Bochum). Support 
for ESS comes from NASA grant NNX09AC14G and for TS comes from DOE grant UD-FG02- 
91ER40626. 



REFERENCES 

Abbasi, R. U., et al. (HiRes-Coll.), Phys. Rev. Letters 104, id. 161101 (2010); arXiv:0910.4T84l 
with a correction in Phys. Rev. Letters 104, id. 199902 (2010). 

Abbott, D. C, Bieging, J. H., k Churchwell, E., Astrophys. J. 280, 671 - 678 (1984) 

Abraham, J., et al. (Auger-Coll.), Phys. Rev. Letters 104 id. 091101 (2010): larXiv:1002.0699l 

Adriani, O., et al. (Pamela Coll.), Nature 458, 607 - 609 (2009); arXiv 0810.4995 

Aharonian, F., et al., (H.E.S.S.-Coll.), Astron. & Astroph. 508, 561 - 564 (2009): llarXiv:0905.0105l 

Ahn, H.S. et al. (CREAM-ColL), Astrophys. J. 707, 593 - 603 (2009): larXiv:091 1.1889 1 

Ahn, H.S. et al. (CREAM-ColL), Astrophys. J. Letters 714, L89 - L93 (2010): la7XTv:1004.1123l 

Alcaraz, J., (AMS-Coll.) et al., Phys. Lett.B 494, 193 - 202 (2000) 

Barvainis, R., Mcintosh, C, Predmore, C. R., Nature 329, 613 - 615 (1987) 

Biermann, L.F., Zeitschr. fiir Astrophys. 29, 274 (1951) 



Biermann, PL., Astron. & Astroph. 271, 649 (1993) - paper CR-I; astro-ph/9301008 



Biermann, P.L., 23rd ICRC, in Proc. "Invited, Rapporteur and Highlight papers"; Eds. D. A. Leahy 
et al., World Scientific, Singapore, p. 45 (1994) 

Biermann, P.L., Proc. Nuclear Astrophysics meeting at Hirschegg, in Proc, GSI, Darmstadt, Ger- 
many, p. 211 - 222 (1998) 

Biermann, PL., et al., Astron. & Astroph. 369, 269 - 277 (2001) 

Biermann, P.L., Bisnovatyi-Kogan, G., & Moiseenko, S., Proc. Brazil Nov 2004 meeting, Magnetic 
fields in the Universe: From Laboratory and Stars to Primordial Structures, Ed. E. d. Pino 
Gouveia et al, AIP Proc. 784, 385 - 395 (2005) 



Biermann, P.L., et al., in "High- Energy Gamma-rays and Neutrinos from Extra-Galactic Sources", 
Heidelberg Jan 2009, Int. J. of Mod. Phys. D 18, 1577-1581 ( 2009a): larXiv:0904. 15071 

Biermann, P. L., Becker, J. K., Meli, A., Rhode, W., Seo, E.-S., & Stanev, T., Phys. Rev. Letters 
103, 061101 f2009b): larXiv:0903.4048l 

Biermann, P.L., Becker, J.K., Caceres, G., Meli, A., Seo, E.-S., & Stanev, T., Astrophys. J. Letters 
710, L53 - L57 (2010); larXiv:0910.lT97l 

Bisnovatyi-Kogan, G. S., Astron. Zh. 47, 813 (1970) 

Bisnovatyi-Kogan, G. S., Moiseenko, S. G., Chinese J. of Astron. & Astroph. Suppl. 8, 330 - 340 
(2008) 

Chang, J., (ATIC-Coll.) et al. Nature 456, 362 (2008) 

Churchwell, E., Bieging, J. H., et al., Astrophys. J. 393, 329 - 340 (1992) 

Dobler, G., Finkbeiner, D. P., Cholis, I., Slatyer, T. R., Weiner, N., Astrophys. J. 717, 825 - 842 
(2010); larXiv:0910.4583l 

Drake, St. A., et al. Astrophys. J. 322, 902 - 908 (1987) 

Ellison, D.C., & Double, G.P., Astropart. Phys. 22, 323 (2004) 

Finkbeiner, D.P., Astrophys. J. 614, 186 - 193 (2004); arXiv:0311547 

Gallant, Y.-A., & Achterberg, A., Month. Not. Roy. Astr. Soc. 305, L6 - L10 (1999) 

Gopal-Krishna, Biermann, P.L., de Souza, V., Wiita, P.J., Astrophys. J. Letters 720, L155 - L158 
(2010); larXiv: 1006.50221 

Heger, A., et al., Astrophys. J. 591, 288 - 300 (2003) 

Hooper, D., et al., Phys. Rev. D 76, ms. 083012 (2007) 



-14- 

Jokipii, J. R., Astrophys. J. 313, 842 - 846 (1987) 

Kardashev, N.S., Astron. Zh. 39, 393 (1962); translated in Sov. Astron. 6, 317 (1962) 

Kardashev, N. S., Astron. Zh. 41, 807 (1964); translated in Sov. Astron. 8, 643 (1965) 

Kolmogorov, A.N., Dokl. Akad. Nauk SSSR 30, 299 - 303 (1941a); 31, 538 - 541 (1941b); 32, 19 
21 (1941c) 

Meli, A. & Quenby, J., Astropart. Phys. 19, 637 (2003a); 19, 649 (2003b) 



Meli, A., & Biermann, P. L., Astron. & Astroph. 454, p. 687 - 694 (2006); astro-ph/0602308 



Moiseenko, S. G., Bisnovatyi-Kogan, G. S., Ardeljan, N., in Astronomy and Beyond: Astrophysics, 
Cosmology, Radio Astronomy, High Energy Physics and Astrobiology. AIP Conf. Proc. 1206, 
282 - 292 (2010) 

Parker, E.N., Astrophys. J. 128 664 (1958) 

Prantzos, N., Adv. Space Res. 4, 109 - 114 (1984) 

Prantzos, N., in Proc. Wolf-Ray et Stars and Interrelations with Other Massive Stars in Galaxies:, 
Eds. Karel A. van der Hucht and Bambang Hidayat, Kluwer, Dordrecht. I.A.U. Symp.143, 
550 (1991) 

Prantzos, N, , in Proc. IAU Symposium 268 Light Elements in the Universe, Eds. C. Charbonnel, 
M. Tosi, F. Primas and C. Chiappini, I.A.U. Symp. 268, in press ( 20101: !arXiv:1003.2317l 

Rickett, B. J., Annual Rev. of Astron. & Astrophys. 15, 479 - 504 (1977) 



Seemann, H., k Biermann, P.L., Astron. & Astroph. 327, 273 (1997); astro-ph/9706117 



Seo, E.S., et al. (CREAM-coll.), Adv. Space Res. 42, 1656 - 1663 (2008) 



-15- 

Stanev, T., Biermann, P.L. & Gaisser, T.K., Astron. & Astroph. 274, 902 (1993) - paper CR-IV; 
astro-ph/9303006l 

Su, M., Slatyer, T. R., Finkbeiner, D. P., submitted (2010): larXTv: 1005.54801 
Volk, H.J., Biermann, P.L., Astrophys. J. Letters 333, L65 - L68 (1988) 
Woosley, S. E., Heger, A., Weaver, T. A., Rev. Mod. Phys. 74, 1015 - 1071 (2002) 



This preprint was prepared with the AAS IATgX macros v5.2. 



